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5 Optical Imaging and Oximetry of Tissue 

CLAIM OF PRIORITY 

This application claims priority mider 35 USC §119(e) to U.S. Provisional Patent 
Application Serial No. 60/364,239, filed on March 13, 2002, the entire contents of which 
is hereby incorporated by reference in its entirety. 

10 TECHNICAL FIELD 

TMs invention relates to optical methods, and more particularly to optical 
metliods of imaging and oximetry. 

BACKGROUND 

Optical techniques provide a safe and non-invasive means of studying human 

1 5 tissue. For example, optical mammography is a technique for detecting cancer in the 
human breast. This technique uses wavelengths of light in the near infrared (NIR) 
spectrum to provide an image where contrast is based on the absorption of the light by 
hemoglobin. Hemodynamic changes such as an increase in the vascular density 
(resulting from angiogenesis) and chaiiges in the blood flow and oxygenation typically 

20 occur at timior sites. See, e.g., S. Zhou et al, Proc. SPIE 2979 :98-106, 1997. The 
increase in vascular density induces a local increase in the optical absorption that may 
allow for the detection of breast lesions, wliich can be benign or malignant. 

The oxygenation level of a lesion has been said to indicate whether a lesion is 
benign or malignant. Measui-ements of the partial pressure of oxygen in tumors have 

25 shown that hypoxic or anoxic conditions often exist in malignant tumors, but not in 

benign lesions. See, P. Vaupel et al., Cancer Res. 5i, 3316-3322, 1991; P. Hohenberger 
et al., Breast Cancer Research and Treatment 48, 97-106, 1998; and P. Vaupel et al., 
Seminars in Oncology 28, 29-35, 2001. Therefore, measurements of the oxygenation 
levels of tissues in vzvo can distinguish benign from malignant lesions in the human 

30 breast. 
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5 SUMMARY 

The invention features improved systems and tecliniques for imaging a region of 
tissue using optical methods and improved techniques for measuring the oxygenation 
level of regions of interest in tissue in a living subject. These tecliniques are non- 
invasive, using safe levels of optical radiation. The optical imaging systems and 

1 0 tecliniques can be applied to a wide variety of samples, including human brain, breast, 
and muscle tissue, as well as any organ that is optically accessible. These techniques 
provide enlianced image contrast, facilitating detection of lesions and tumors and 
differentiating between benign and malignant varieties. In addition, the new tecliniques 
for measuring the oxygenation level of a region of interest, such as a lesion or tumor, are 

15 useful, e.g., to further assess the malignancy of these regions. 

hi general, the new systems and techniques improve the diagnostic capability of 
optical maimnography. While these systems and methods can be used to image breast 
tissue and oximetry of breast tumors, they can also be used to measure the oxygenation 
level of other hemoglobin-rich localized tissue areas. For example, the focal increase in 

20 cerebral hemoglobin concenti'ation induced by selected cerebral activity, a localized 

hematoma, and large blood vessels can be investigated using the methods presented here 
to measure oxygenation levels. 

In one aspect, the invention features imaging systems including one or more light 
sources that provide light to a sample, a detector that detects light transmitted through the 

25 sample, and a processor that calculates data using spatial second derivatives of products 
of thiclaiesses of the sample and mtensities of the transmitted Ught. The new systems can 
include a number of features. For example, the light sources can be laser diodes, the 
detectors can be photomultiplier tubes or avalanche photodiodes, and the processor can 
be a computer. The systems can be implemented using knovra equipment, and the light 

30 sources and detectors can be separated by a fixed distance, e.g., to simpHfy data analysis. 
In addition, the systems can include a monitor to display an image of the data and/or a 
computer-readable medium to store the data. 
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5 The invention also features a computer-readable medium having a program that is 

used by a processor to receive data from a sample that represents measures of thiclaiesses 
of the sample and intensities of light transmitted tlirough tlie sample, and to process the 
data from the sample using spatial second derivatives of products of tlie tliicknesses of 
the sample and intensities of the transmitted liglit to generate values tliat represent an 

1 0 image of the sample. This computer-readable medium is useful to transport and 
distribute the new algorithms to various users. 

In a further aspect, the invention features a computer-readable medium including 
a progi-am that is used by a processor to receive data from a sample that represent 
measures of intensities of light transmitted through the sample for at least two 

1 5 wavelengths of light, select a pafr of wavelengths to minimize a difference betv,'een 

relative changes in the intensities of the light transmitted through a region of the sample, 
where the relative changes in intensities of the light are measured relative to a 
backgi-ound intensity of light transmitted through the sample, and process the relative 
changes in intensities of the light at the pair of wavelengths using measures of a 

20 background reduced scattering coefficient of tlie sample at the pair of v^favelengths and 
molar extinction coefficients for deoxy-hemoglobin and oxy-hemoglobin to generate 
values that represent an absolute level of oxygenation of the region of the sample. 

The computer-readable medium can also include both programs to process spatial 
second derivatives and to calculate absolute levels of oxygenation of a sample, as 

25 described herein. 

In other aspects, the invention features metliods that can be carried out using the 
new systems, such as methods of detecting at least one region of a sample having an 
absorption level different from a background level of absorption in the sample. These 
methods involve obtaining thiclaiesses of the sample and intensities of light transmitted 

30 thi-ough the sample at a plurality of locations, calculating data points for the locations 
using spatial second derivatives of products of the thiclaiesses of the sample and the 
intensities of the transmitted light, and comparing the data points to detect a region of the 
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5 sample having an absorption level different from a backgromid level of absoiption within 
the sample. 

Using the spatial second derivatives enhances differentiation of regions of the 
sample having different levels of absorption from the background level of absoiption 
within the sample, hi addition, this approach has the advantages of being insensitive to 

1 0 the size, shape, and depth of the region within the sample. Embodiments of tliis aspect of 
the invention include one or more of the following features. The sample can be tissue, 
e.g., breast, brain, or muscle tissue. Consequently, the methods and systems are useful in 
performing optical mammography, studying brain activation, and identifying ischemic or 
underperfused regions in muscle. The region of the sample having an absorption level 

15 different from the background level of absorption witliin the sample can be a lesion or 
tumor. This region can be cancerous. Thus, the systems and methods can be applied to 
identify lesions, including cancerous tumors, witliin breast tissue. The systems and 
methods can also be used to identify tumors regardless of their size, shape, or depth of the 
tumor within the breast. 

20 The methods can fiirther involve illuminating the sample with Ught and detecting 

light transmitted through the sample. This light can have two or more wavelengths. 
Illumination with multiple wavelengths of light is usefid to apply the new methods of 
oximetry discussed herein. One or more of these wavelengflis can be in the near infrared 
region of the spectrum. Infrared light exhibits a desirable amount of absorption within 

25 human tissue, e.g., breast tissue. The transmitted light can be generated by coplanar 
sources and the plurality of locations where the light is detected can be coplanar. Using 
coplanar sources and detectors facilitates data analysis. 

Li these systems and methods, the data points N"{x,y) are calculated using the 
second derivative of an equation N{x,y) - roaco/r{x,y)ac{x,y), where x and ai'e respective 

30 X- and y-coordinates of the coplanar locations, rg is the maximum thickness of the 

sample, aco is the intensity at the location where the thicbiess of the sample is ro, r{x,y) is 
the thicloiess of the sample at location {x,y), and ac{x,y) is the intensity at location {x,y). 
These data points can be displayed in an image, as this fomiula provides useftil data 
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points for imaging. Alternatively or additionally, they can be stored in a computer- 
readable medium. By storing these data points on a computer-readable medium, the 
image can be used for later comparison or displayed at alternate locations and times. 

hi another aspect of the invention, method are used, in the new systems, to 
calculate an absolute oxygenation level of a region of a sample by selecting a pair of 
wavelengths of light to minimize a difference between relative changes in intensity of the 
light transmitted tlurough the region of the sample, where the relative changes in intensity 
of the light are measured relative to a backgi-ound intensity of light ti-ansmitted through 
the sample; obtaining measures of background reduced scattering coefficients of the 
sample at the pair of wavelengths; and calculating tlie absolute oxygenation level of tlie 
region of the sample using the relative changes in intensity of the light at the pair of 
wavelengths, the measures of the background reduced scattering coefficients of tlie 
sample at the pair of wavelengtlis, and molar extinction coefficients for deoxy- 
hemoglobin and oxy-hemoglobin. 

These methods provide accurate quantitative measurements of the oxygenation 
level of regions of interest witliin a sample. Such measurements are valuable to 
differentiate between benign and cancerous tumors within breast tissue. Using an 
appropriate pair of wavelengths of hght provides measurements that ai-e largely 
independent of the tumor size, shape, and location inside the breast. It is significant that 
the pah of wavelengths for the near-infrared measurement of tumor oxygenation is not 
fixed, but instead depends generally on tumor oxygenation. 

These methods (and the computer readable medium) can include the use of the 
formula 

^ 1^.0(^2 ) M/Jf^^^^ 



[SHb(^2)-£Hb02(^2)] + [£Hb02(^l)-eHbai)]4^^^^^y^ 
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5 where SO2 is the oxygenation level of the region, M I ||^^^ is a maximum value of the 
relative change in intensity at a fu-st wavelength in the pair of wavelengths, hi 1 1^"-^ is 
a maximum value of the relative change in intensity at a second wavelengtli in the pair of 
wavelengths; - ^f is a ratio of the background reduced scattering coefficient at the 

first wavelength to the background reduced scattering coefficient at the second 

1 0 wavelength, 8Hb(^i) and 8Hb(^2) are the molar extinction coefficients for deoxy- 

hemoglobin at the first and second wavelengths, and SHb02(?^i) and ZmoiQ^i) are the molar 
extinction coefficients for oxy-hemoglobin at tlie first and second wavelengths. Tliis 
fonnula provides accurate values of oxygenation levels and is useful for ascertaining the 
malignancy of tumors. The difference between the relative changes in intensity can be 

15 zero, as this provides an optimal pair of wavelengths. 

These methods can also include a number of additional features. For example, the 
sample can be illuminated with two or more wavelengths of light and the light 
transmitted tluough the sample can be detected. The method can also include displayhig 
the absolute oxygenation level of the region in an image and/or storing the absolute 

20 oxygenation level of the region in a computer-readable medium. Storage of oxygenation 
level infonnation facilitates transport and sharing of this information, as well as 
permitting analysis at another location and/or time. 

In another aspect, methods are used to determine whether a tumor in a tissue 
sample is a malignant tumor or a benign tumor. Tlnclmesses of the sample and mtensities 

25 of light transmitted tlirough the sample at a plurality of locations for two wavelengths of 
light are obtained. Spatial second derivatives of products of the sample thiclaiesses and 
the intensities of the ti-ansmitted light at the locations for the two wavelengths of light are 
calculated. An oxygenation level of the tumor is then calculated based on the spatial 
second derivatives for the two wavelengths of Hght, the molar extinction coefficients of 

30 oxy-hemoglobin for the two wavelengths of light, the molai' extinction coefficients of 
hemoglobin for the two wavelengths of light, relative changes in intensity of the light for 
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5 the two wavelengths and the measures of the background reduced scattering coefficients 
of the sample for the t^vo wavelengths. The oxygenation level of the tmnor is compared 
with the oxygenation level of non-tumor regions of the sample to deteraiine whether the 
tumor is malignant or benign. 

These methods provide accurate measurements of the relative oxygenation level 
1 0 of regions of interest witliin a sample. Such measurements are valuable to differentiate 
between benign and cancerous tmnors within, for example, breast tissue, brain tissue, and 
muscle tissue. 

These methods may use the following formula for calculating an oxygenation 

level: 

15 

^^^_^[HbOJ_, where 
A[HbOJ'+A[Hb]* 



[^2^"(X,)8„b02(^/)j 




hi 


T N'Ui)^m SHb02 (:^/)SHb (.h)] 


1 
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2EHb02(^/)EHb(^/)j 


2 



A[Hb]* = 







i-l 




j^ZeHb02(^/)EHb(^/)l 


1 


'l:^Hb02(^/)] 




hi 




Hb(^,)~^ 


2 

1 



and where OL is the oxygenation level of the tumor; i is a wavelength index for the two 
wavelengtlis; SHb02 and Snb are the molar extinction coefficients of oxy-hemoglobin and 
deoxy-hemoglobin, respectively; A[Hb02]* and A[Hb]* are relative values for tlie spatial 
changes in the concentrations of oxy-hemoglobin and deoxy-hemoglobin, respectively; 
25 and 

N" is a spatial second derivative of an intensity of transmitted light. 
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5 These methods can include one or more of the following features. The sample 

may be illuminated with a plurality of wavelengths of light, and the light transmitted 
tkough the sample may be detected at a plurality of locations. An image of oxygenation 
levels at the plurality of locations within the tissue sample may be displayed. Values of 
oxygenation levels at the plurality of locations within the tissue sample may be stored in 

1 0 a computer-readable medium. 

In another aspect, the invention features methods for diagnosing cancer in a 
sample by obtaming thicknesses of tlie sample and intensities of light transmitted through 
tlie sample at two or more of locations; calculating data points for tlie locations using 
spatial second derivatives of products of the thiclaiess of the sample and the intensities of 

15 the transmitted hght; comparing the data pomts to detect a region of the sample having an 
absoxption level different from the background level of absorption within the sample; 
selecting a pair of wavelengths of light to minimize the difference between relative 
changes m intensitj' of the light transmitted tlirough the region of the sample, where the 
relative changes in intensity of the light are measured relative to the background intensity 

20 of light transmitted thi-ougli the sample; obtaining measures of background reduced 
scattering coefficients of the sample at tlie pair of wavelengths; calculating an absolute 
oxygenation level of the region of the sample using the relative changes in intensity of 
the hght at the pair of wavelengths and tlie measures of the background reduced 
scattering coefficients of the sample at the pair of wavelengths; and diagnosing whether 

25 the region of the sample contains cancer based on the absolute oxygenation level of the 
region. 

Unless otherwise defined, all technical and scientific terms used herein have the 
same meaning as commonly understood by one of ordinary skill in the art to which this 
invention belongs. Although methods and materials similar or equivalent to those 
30 described herein can be used in the practice or testing of the present invention, suitable 
methods and materials are described below. All pubhcations, patent applications, 
patents, and otlier references mentioned herein are iacorporated by reference in their 
entirety. In case of conflict, the present specification, including definitions, will control. 

-8- 
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5 111 addition, the materials, methods, and examples are illustrative only and not intended to 
be limiting. 

The details of one or more embodiments of the invention are set forth in the 
accompanying drawmgs and the description below. Other features, objects, and 
advantages of the invention will be apparent from the description and drawings, and from 
10 the claims. 

DESCRIPTION OF DRAWINGS 

FIG. 1 is a schematic diagram of a system of the present invention. 
FIG. 2 is a prototype of the new system used to obtain optical data for second- 
derivative images. 

15 FIG. 3 is a prior art image of the left breast, craniocaudal projection, for a 58-year 

old patient affected by invasive ductal carcinoma 

FIG. 4 is a line graph derived from the data shown in FIG. 3 at>' = 3.2 cm. 

FIG. 5 is second-derivative image of the data generated by the new methods. 
FIG. 6 is a line graph derived fi'om the data shown in FIG. 5 at;-- = 3.2 cm. 
20 FIG. 7 is a schematic diagram of a breast shghtly compressed between two 

parallel plates. 

FIG. 8 is a schematic diagram of an optical arrangement for evaluating optical 
properties in test samples. 

FIG. 9 is a diagram of an apparatus for measuring optical signals in test samples. 
25 FIG 1 0 is a table of the results of the measurements of the optical properties of 

inclusions and background medium in a test sample. 

FIG. 11 is a graph of optical intensity (I) measured during a linear scan across a 
test sample. 

FIG. 12 is a graph of experimental results for - ^ / -^olmax ^ function of the 
30 product fj-io' A(ifl for irregularly shaped inclusions m a test sample. 

FIG. 13 is a graph of experimental results for - A/ /J'oLax ^ function of the 
product |J.io'A|i„ for a cylindrical inclusion in a test sample. 
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FIG. 14 is the calculated dependence of - A/ //o^^^^ on the product |a,o'A|i„ for a 
spherical object embedded m a miifonn turbid medium. 

FIG. 15 is a gi-aph showing the dependence of A/ / 1^]^^^ on fi„o for a test sample. 

FIG. 16 is a spectrum of the product |iio'A|a« for SO2 values of 0, 20, 40, 60, 80, 
and 100%. 

FIG. 17 is graph obtained using prior art perturbation analysis for a number of 
wavelength pairs for a 1.4 cm-diameter sphere (on the mid-line between som-ce and 
detector) having a range of oxygenation levels. 

FIG. IS is graph obtained using prior art perturbation analysis for a number of 
wavelength pairs for a 3.0 cm-diameter sphere (on the mid-line between source and 
detector) having a range of oxygenation levels. 

FIG. 19 is graph obtained using prior art perturbation analysis for a number of 
wavelength pairs for a 1.4 cm diameter sphere (off the mid-Hne by 1.5 cm) having a 
range of oxygenation levels. 

FIG. 20 is a graph obtained using the new methods for a selected wavelength pair 
for a 1 .4 cm-diameter sphere (on the mid-line between source and detector) having a 
range of oxygenation levels. 

FIG. 21 is a graph obtained using the new methods for a selected wavelengtli pair 
for a 3.0 cm-dianieter sphere (on the mid-line between source and detector) having a 
range of oxygenation levels. 

FIG. 22 is a graph obtained using the new methods for a selected wavelengtli pair 
for a 1.4 cm diameter sphere (off the mid-line by 1.5 cm) havuig a range of oxygenation 
levels. 

FIG. 23 is a graph of S02^'^ versus A[x„{Xi) for SO2 values of 0, 20, 40, 60, 70, 80, 
and 100%. 

FIGs 24al to 24c2 are a series of images of a ductal carcinoma in a left breast. 
FIGs. 25a to 25c are a series of images of a benign mastopathy in a left breast. 
FIGs. 26a to 26c are a series of images of benign inflammation in a breast. 
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5 Like reference symbols in the various drawings indicate like elements. 

DETAILED DESCRIPTION 
The invention features optical systems and methods for imagmg based on 
oxygenation levels and quantifying oxygenation levels within a region of interest. These 
systems and methods are useful to assess differences in oxygenation levels within 

1 0 mammalian, human, e.g., tissue and have particular application to identifying cancerous 
tumors in breast tissue in vivo. The imaging techniques provide enhanced contrast by 
mapping the second derivative of the optically acquired data. Measurement of absolute 
oxygenation levels involves selection of an appropriate pair of wavelengths to minimize 
the difference between the relative change in the intensity of the light detected at those 

15 wavelengths. 

Optical Systems 

An optical system for scanning a human breast according to the present invention 
is shown in FIG. L Two glass plates 10 are used to slightly compress the breast. Two 

20 distance sensors 20 and 30 positioned next to the illumination 40 and collection fibers 50 
measure the free space between tlie compression plates and the breast. The difference 
betv^'een the plate separation and this free space yields the breast thickness at each pixel. 
The distance between the two plates is adjustable to fit different breast sizes, and it is 
possible to rotate tlie whole assembly to allow for the acquisition of breast images in the 

25 craniocaudal, oblique, and mediolateral projections. The output of the distance sensors 
(or, as an alternative, the optical signal detected by the scanning collection fiber 50) is 
also used to detect the proximity to the breast edge. By driving the scanner to the next 
line using the driving circuit 60 when the detected breast thickness becomes less than 1 
cm (or when the optical transillumination signal tlu'ough the breast exceeds a given 

30 threshold), the system only scans the area actually occupied by the breast. 

The light source can be a Xe arc lamp 70 (ORIEL histruments. Model No. 6258), 
and the light is passed tlirough a band-pass filter 80 to emit over the specfral range 680- 
880 mn. This continuous wave (CW) light source provides a complete range for 
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5 selection of optimal wavelengths. The lamp 80 is coupled through a F/1 .0 condenser and 
a converging lens (collectively shown as 90) to either or both of the 3-5 mm diameter 
glass fiber bundles 40 and 100. A fiber bundle 100, which is a fixed illumination fiber, 
can be used for background spectroscopy. The collection fiber 50, a fiber bundle with an 
internal diameter of 3-5 nun and a numerical aperture of 0.22, delivers the light to the 
1 0 entrance slit of a 1 5 cm-f74 spectrogi-aph 1 1 0 (Roper Scientific, VersArray: 5 1 2B) for the 
parallel measurement of the whole spectnim at each pixel. The CCD camera chip 120 is 
appropriately binned to achieve the best compromise between signal-to-noise ratio and 
spectral resolution. The spectral data for the whole image is stored in a computer 130 
and is post-processed. A display 140 displays an image of the post-processed data. 

15 

Second-Derivative hnaging 

Processing optical data from breast tissue using spatial second derivatives further 
enliances the detection of areas of higher light absorption m the breast. Light absorption 
is a function of concentration and oxygenation levels. The second-derivative images at 

20 multiple wavelengths can be used in spectral analysis, which improves the specificity of 
optical mammography by discriminating between benign and malignant breast lesions. 

Optical data for second-derivative images were acquired using a prototype 
designed and built by Siemens Medical Engineering (Erlangen, Germany) shown in FIG 
2. Other devices can be used to acquire images. This prototype did not feature a CW 

25 light source, histead, four laser diodes (690, 750, 788, and 856 nm) modulated at a 

fi-equency of -70 MHz were used as optical sources. The light fi-om these diodes was fed 
into optical fibers, wliich were joined to a fiber coupler. The source and detector fibers 
are located on opposite sides of the breast, which was slightly compressed between two 
parallel glass plates. The light received from the detector fibers passed through a 

30 photomultiplier tube and analyzer to a computer. Amplitude and phase images were 
obtained by scamiing the source and detector fibers in tandem. Two projections of each 
breast were typically acquired, craniocaudal and oblique. The raw data was con-ected for 
edge effects arising firom tiie variable thickness of the breast between the two plates. The 
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5 edge-coiTected data yield was used to create images based on an N-parameter, defined as 
N(x,y) = r,acjiix,y)ac{x,y), -where tq is the maximum thiclmess of the sample, aco is the AC 
amplitude at a pixel where the breast thickness is ro, acix,y) is the amplitude measui-ed at 
pixel {x,y), and r{x,y) is the breast tliiclcness at that pixel derived from the phase 
information, see S. Faiitini et al., Med. Phys. 23:149-57, 1996. A representative A^-image 

10 of the left breast, craniocaudal projection, for a 58-year old patient affected by invasive 
ductal carcinoma is reported in FIG. 3 . A line graph ofN <Ay=3>.2 cm is shown in FIG. 
4. 

FIG. 5 is an image obtained by taking the second derivative {N"). Tliis image was 
generated by first smootliing the A^-image and then taking the sum of the second 
15 derivatives at each pixel {ij) calculated in the vertical (y), horizontal (;c), and two (x+v 
and x-y) diagonal directions, which are defined as follows: 



(A,)^+(A,)^ 



(A,.)^+(A,.)^ 

where i and j are the pixel indices in the and j directions, respectively, while A,- and A,- 
are the pixel size in the x and y directions, respectively. Altematively, instead of the sum 
one can take the maximum value of the second derivative. 

FIG. 6 shows a line graph ofN" at y = 3.2 cm. In FIG. 5, areas with positive TV" are 
set to white (by setting a thi-eshold at N" = 0), areas with negative N" (corresponding to 
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5 regions of high optical absorption, and generally indicative of lesions) are displayed in 
gi-ay-scale, and the region external to the breast is set to black. This image has enhanced 
contrast and identifies those areas in the image characterized by a local maximum. As 
FIG. 5 demonstrates, the secondary peak at x = 9 cm is now greatly enlianced, and an 
even smaller feature around x = 6 cm is now visible. 

10 

Measurement of Oxygenation Levels 

To a first approximation, a locaUzed optical perturbation embedded in a tm-bid 
medium can be modeled using the fu-st-order perturbative solution to the diffusion 
equation. First-order perturbation appUes to small optical perturbations, where "small" 

15 means that the linear dimensions of the perturbation ai-e much smaller than the distance r 
between the illumination and the collection points that induce a small (with respect to 1) 
relative change in the optical signal. Tliis model is schematically illustrated in FIG. 7 for 
a lesion embedded in a transilluminated human breast, where the lesion is located at a 
distance ro fi'om a glass plate on the detector side and at a distance r-j'o fi-om a glass plate 

20 on the source side of the breast model. 

The absorption and reduced scattering coefficients of the backgi-ound medium 
(healthy breast tissue) are denoted with |j,„o and \Xso', respectively, and the absorption and 
scattering perturbations (associated with the tumor) with A\x„ and A|i;, respectively. This 
means that the absoiption and reduced scattering coefficients at the location of the optical 

25 peiturbation (tumor) are given by n„o + Aii„ and |j,^o' + A^/, respectively. The relative 
intensity change AI/Iq (where /q is the unperturbed intensity measured in the background 
medium) induced by the optical perturbation is maximized (in absolute value) when the 
source, detector, and the perturbation are coUinear as shown in FIG. 7. The first-order 
perturbative solution to the diffusion equation leads to the following expression for this 

30 maximal intensity change, see Boas et al., Appl. Opt. 36:75-92, 1997, which is 
incorporated by reference herein: 
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(1) 

V is the volume of the lesion, and the geometrical parameters r and rg are defined 
in FIG. 7. The second term inside the square bracket in Eq. (1), the scattering 
pertuitation temi, is neglected on the basis of in vivo studies tliat have shown that the 
1 0 scattering tumor-to-background contrast is much smaller than the absoiption contrast, 
see, e.g., Fantini et al, App. Opt. 37:1982-89, 1998, which is incorporated herein by 
reference. Therefore, Eq. (1) reduces to: 




As shown by Eq. (2), first-order perturbation theory predicts that the maximal 
1 5 intensity effect of an optical inclusion (for the case An/ = 0) is given by the product of a 
wavelength-independent geometric factor and a wavelength-dependent factor given by 
|Li.o'A|.i„. 

hi the new oximetry techniques, Eq. (2) is generalized to predict the functional 
dependence of A7 /^oLax ^V^" beyond tlie Umits of vahdity of the first- 

20 order peiturbation analysis. Specifically, in the case Afi/ = 0, Eq. (2) is generalized to 
the case of spatially extended objects that induce significant changes in the optical 
intensity. This is accomphshed by expressing AI //oLax ^ fiinction of the 

product \iso'A[ia and wavelength-independent parameters that are related to the source- 
detector separation (?-), and to the object-detector distance (ro), size (F), and shape (s) of 

25 tlie object. Therefore, AI //olmax expressed as: 




= f{r,lo,V,S,[\^'sA)^a]) 



(3) 
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5 where / indicates an unknown function of its arguments. 

As in the perturbation case, AI I Iq \^^^ is a monotonic function of the product 
|a.io'A|an. The wavelength dependence of A/ //oLax ^^'^^^ appears impHcitly in the 
ai-gument \Xso'^\in, because the geometrical pai-ameters 7-0, V, and s are mdependent of 
the wavelength. Consequently, because of the monotonicity of A/ / /q [^^^ on |j.jo' Aia„, if 

1 0 two wavelengtlis X\ and X2 are selected such that AI I I^f^^. =AI/Iq f^^^ , then 

\i.so'ih)^\i-a{h) = )aio'(^2)Aix„(A,2). Therefore, the ratio of the absorption perturbations at 
these two wavelengths is given by the inverse of the ratio of the background reduced 
scattering coefficients at the same two wavelengths: 

(4) 

This analysis indicates that by appropriately choosing the two wavelengths X] and 
X,2, one can translate a measurement of the background scattering ratio into a 
measurement of the absorption-perturbation ratio associated with the embedded 

20 inclusion. This result is relevant for tlie near-infirared determination of the oxygen 
saturation of hemoglobin in breast lesions because (1) the near-infirared spectral 
properties of hemoglobin and breast tissue are such that it is usually possible to identify 
two wavelengths at which the tumor-induced intensity perturbations are sindlar, and 
because (2) the oxygen saturation of hemoglobin is only a function of the ratio of the 

25 optical absorption at two wavelengths, see Fantini et al, Opt. Eng., 34:32-42, 1995. If 
the absorption at the tumor location is given by the sum of the contributions from the 
background {\i„o) and fi-om the tumor (|j.„^''), then Ap,„ can be identified with the tumor 
absorption {\xj-^^ = A|a,„). If, instead, the absoiption at the tumor location is only due to the 
tumor, then n„^'^ = \iao + A(x„. However, if A^a » |a„o (high-contrast tumor), it can still be 

30 maintained that A|j,„ is representative of the tumor absorption » An„). 
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Here, the oxygen saturation of hemoglobin (SO2) associated witli the additional 
absorption A[J,„ is considered. At least in the two cases mentioned above, SO2 is 
representative of the tumor oxygenation. By assuming that A\x„ is solely due to 
hemoglobin, the expression for SO2 in tenuis of A|a„(Xi) and A^„(A,2) becomes the 
following: 



where SHb and EHb02 are the molar extinction coefficients of deoxy-hemoglobin and oxy- 
hemoglobin, respectively. 

In practice, it may not be possible to identify two wavelengths A-i and X2 that 
exactly satisfy the requirement that A7 / /q |J^^J = AI I I^f^J^ . In this case, the optimal 
pair of wavelengths (X,i,A,2) minimizes the absolute value of the difference 
MIiJl'I - AI /Iofi£ . To coixect, at least in part, for a difference between AZ/AP'^ 

I max I max ui^ax 

and AI //olll'^ ' dependence AI /iolmi °^ l^io('*^)Al^n(^) given by Eq. (2) is used to 



SO2 = 



gHb (^2) - SHb (^1 )A^, ^2)/^}!, (h) 



[e Hb (>^2 ) - £ Hb02 (>^2 )] + [e Hb02 (^1 ) - £ Hb C^i )] AH„ ^2 )/A^ „ (5^1 ) ' 



write: 



A^„(?^i) ^'.0(^2) Al/lf 

In 




I max 
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5 By combining Eqs. (5) and (6), the lesion saturation is given by: 



(7) 

which is the equation used to quantify the oxygenation level of, e.g., a tumor, from 
1 0 measurements of (1) the background (healthy tissue) reduced scattering coefficient 
[M'jo'WIj and (2) the maximal relative change in the intensity caused by the tumor 

(A///o£l)- 

This approach is not exclusively based on first order perturbation theory because 
the criterion to select the two wavelengths and X,2 refers to situations that are beyond 
1 5 the limits of appUcability of first-order perturbation theory. As a result, even if 

A///o|j^g^ does not show a linear dependence on |J.jo'Afx„ as predicted by first-order 
perturbation theory, these methods still provide accurate readings of oxygenation levels, 
provided that two wavelengths A,i and ^2 can be identified such that 

20 In addition to seelcing a pair of wavelengths that minimize the intensity change 

induced by the mclusion, it is preferable to use wavelengths that are not too close to one 

another. The sensitivity of the oxygenation measurement is degraded if the two near- 
infrared wavelengths are too close to each other, see Sevick et al, Anal. Bioclim. 
195:330-51, 1991. To avoid this degradation, the criterion that |X,2 - >x canbe 
25 adopted to guarantee that the minimization of the absolute value of the difference 

^ ' ^0 Imax -^Ih Imax ^^^^^ appropriate pair of wavelengths. Possible values of x 
include, e.g., any distance between and including 40 and 60 nm. 
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5 This approach has the potential of being more robust in practical clinical 

measurements with respect to full reconstruction schemes. The only information used by 
the new method is (1) the spectrum of the backgi-ound scattering coefficient, n^o'{X), 
(2) the background intensity, kiX), and (3) the maximum intensity change (Ai)n,ax caused 
by tlie tumor. To best apply this method, the spectral measurements should be conducted 

10 at a large number of wavelengtlis in the range of 680 to 880 nm, for example 
continuously over this spectral band. 

The reduced scattering coefficient of the background medium as a function of 
wavelength, |J,so'(^). can be measured by averaging time-resolved measurements at 
several breast locations. Because of the featureless scattering spectrum, measurements at 

15 a few discrete wavelengths Xi can be effectively extrapolated to yield a continuous 

spectrum of \iso'(X), as demonstrated by Bevilacqua et al., Appl. Opt. 39:6498-507, 2000. 
With respect to the measurement of the background intensity /o(X,) in tlie presence of a 
heterogeneous background such as breast tissue, it can be appropriate to consider an 
average background intensity over a specifically selected breast area rather than the 

20 intensity measured at a particular breast location. 

It is also possible to use the transmitted Ught images to provide information 
concerning the relative oxygenation level of the detected optical inhomogeneities. 
Multiple second-derivative images taken at different wavelengths may be used to 
calculate a relative oxygenation index. In particular, pseudo variations in the 

25 concentrations of oxy-hemoglobin (A[Hb02]*) and deoxy-hemoglobui (A[Hb]*) may be 
determined and then combined to calculate a relative oxygenation index according to the 
following equations: 
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A[Hb]* 



SHb02(^,)EHb(^;) 
[l^Hb02(^/)Jl:4b(^/)]-[ZBHb02(^,-)BHb(^,)j 



[^SAr"(X,)8Hb(^,-)j 


[EeHb02(^/)] 


l-l 


7V"(X,.)eHb02 Eh^,02 (X.OEHbC^.oj 


[i:4b02(^/)][2:sL(^.)] 


l-l 


ZsHb02(^,-)eHb(^/)l 


2 



A[Hb02]* +A[Hb]* 



(10) 



Here, i is the wavelength index, and £Hb02 and SHb are the molar extinction 
coefficients of oxy-hemoglobin and deoxy-hemoglobin, respectively. The stars in 
A[Hb02]* and A[Hb]* refer to indicative (or pseudo) values for the spatial changes in the 
concentrations of oxy-hemoglobin and deoxy-hemoglobin associated with the spatial 
changes in the optical absorption coefficient An„. Tlie actual spatial changes A[Hb02] 
and A[Hb] would be obtained by replacing ¥i"(ki) with A[i„{Xi) in Eqs. (8) and (9). At 
least in the first-order perturbation limit, it can be shown that there is a direct 
proportionality between A^' ' and A^.^, which is the basis for using N' ' in Eqs. (8) and (9). 
However, the proportionality between N' ' and A|aa does not hold for the strong 
perturbation on the optical data caused by relatively large breast tumors. For this reason, 
stars in A[Hb02]* and A[Hb]* are used, the ratio of Eq. (10) may be known as a relative 
oxygenation index rather than oxygen saturation of hemoglobm. 
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Diagnosing Cancer 

Tlie imaging and absolute oxygenation level quantification tecliniques described 
herein can be used in various combinations to diagnose the presence of cancer, for 
example, in human breasts. For example, one can use any one or two of the N, N", or 

1 0 oxygenation level images, or one can use all three images of a given tissue sample. 
Combining these tecliniques, and comparing tlieir results, affords enhanced images 
supported by quantitative measurements that improve diagnostic potential. 

Optical mammography of human breasts is perfonned using a system such as that 
shown in FIG. 2. Two projections of each breast are typically acquired, craniocaudal and 

15 oblique, by illuminating the breasts with multiple wavelengths, preferably with at least 
some of the wavelengths in the range of 680 to 880 rnn. Amplitude and phase unages are 
obtained by scanning the source and detector fibers in tandem, and the raw data was 
corrected for edge effects arising from the variable thickness of the breast between the 
two plates. Second-derivative unages such as FIG. 5 are generated for each breast. 

20 These unages display differences in Ught absorption witli enhanced contrast, and 

inclusions with high absorption, such as lesions, can be identified. For example, in FIG. 
5, the feature located at x = 12 cm is a lesion that represents a region of mterest. 

The oxygenation levels in the identified regions of interest are quantitatively 
assessed. To perform the quantitative assessment, a pair of wavelengths is selected to 

25 minimize the intensity change induced by the region of interest. Since it is preferable to 
use wavelengths that are not too close to one another, a boundary condition of, e.g., 40 
nm of separation in the two wavelengths is applied. Using these two wavelengths, Eq. 
(7) is applied to compute the absolute oxygenation level within the region of interest. 
This absolute oxygenation level improves the ability to distinguish betv\'een maUgnant 

30 and benign inclusions. Because measurements of the partial pressure of oxygen in 

tumors have shown that hj'poxic or anoxic conditions often exist in malignant tumors, but 
not in benign lesions, a low oxygenation level of the region of interest is an indicator of 
its maUgnancy. Differences of 5% to 15% m oxygenation levels between normal and 
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5 malignant tissue have been observed, see, M J. Holboke et al., J. of Bionied. Opt. 5, 237- 
247, 2000; Tromberg et al., Neoplasia 2, 26-40, 2000; Tromberg et al., Philos. Tmas. R- 
Soc. London 352, 661-668, 1997. 

The invention is fiirther described in the following examples, which do not limit 
the scope of the invention described in the claims. 

10 EXAMPLES 
Test Samples and Apparatus 

Inclusions were comprised of a mixture of General Electric silicones, model No. 
RTV61 5 (clear) and model No. RTVl 1 (white). The white sihcone was used as the 
scattering material and black India ink was used as the absorber. Two different mixtures 
15 were prepared with different optical properties, (the first with \i„ ~ 0.05-0.06 cm"' and 
~ 9-10 cm'' and the second with ~ 0.11-0.15 cm"' and [is ~ 9-10 cm"' over the 
wavelength range of 752-840 mn considered) to cover a range of A|j,a that is 
representative of the expected absorption contrast of breast lesions in vivo. 

The optical characterization of the inclusions was perfonned using the 
20 an-angement shown in FIG. 8. The silicone mixtures were formed into two cylindrically 
shaped blocks (one to characterize each mixture) of 1 1 cm diameter and 9 cm height. A 
jfrequency-domain, near-infrared spectrometer (OxiplexTS, ISS, Inc., Champaign, IL) 
housed and controlled the laser sources and optical detectors used to obtain the 
measm-ements. The sources were laser diodes at six discrete wavelengtlas (752, 778, 786, 
25 813, 830, and 840 nm) and the optical detector was a photomultiplier tube (Hamamatsu 
Photonics R928). The sources M'ere mtensity-modulated at a frequency of 1 10 MHz, 
electronically multiplexed at a rate of about 10 Hz to time-share the detector, and were 
coupled to 400 |j,m core-diameter optical fibers that were collected into a fiber bundle 
with rectangular cross-section 1 .2 nmi in width and 3 mm in length. The detector was 
30 coupled to anotlier fiber bundle of circular cross-section, 3 mm internal diameter. The 
detector fiber bundle remained fixed and was in contact with the surface of the silicone 
blocks, while the end of the source bundle was brought close to the surface of the blocks 
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5 (< 0.25 mm away). The starting distance between the source and detector fibers was 1 
cm (x-direction). The source fiber was then moved in the negative x-direction defined in 
FIG. 8 (away from the detector fiber) at a rate of 0.65 mm/s. Data was acquired at every 
0.5 s providing a measurement every 325 \xm over a total traveled distance of 2 cm. The 
amplitude (AC), average value (DC), and phase of the detected modulated mtensity were 

1 0 recorded and translated into measm-emeiits of the reduced scattering and absorption 
coefficients by employing a multi-distance method based on the diffusion equation and 
semi-iiifinite boundary conditions, see Fantini et al., J. Opt. Soc. Am. 11:2128, 1994, 
which is incorporated herein by reference. 

The background medium was comprised of a mixture of 1 liter of Liposyn 10% 

15 (Abbott Laboratories, North Chicago, IL) and 8 liters of deionized water and was 
contained in a rectangular vessel with dimensions of 27 cm width, 32 cm length, and 
13 cm height. Again, the ISS Ixequency-domain spectrometer was used to perform the 
measurements. The experimental arrangement was similar to that used in FIG. 8 except 
the ends of the source and fiber bundles were fully immersed in the Liposyn mixture to 

20 simulate an infmite medium. Acquisition times, laser multiplexing rate and source- 
displacement speed (with the detector remaining fixed) were also the same as in the 
inclusion characterization. The AC, DC and phase data were used in the diffusion model 
for light propagation in an infinite, highly scattering medium to determine the 
background reduced scattering (^.^0') and absorption coefficients (p,„o). Measurements of 

25 |i„o over the wavelengths used ranged fi-om 0.02 cm'' to 0.04 cm"' and those of |j,jo' 
ranged from 10 cm"' to 12 cm"' and were chosen to match typical absorption and 
scattering coefficients in breast tissue for this wavelength range [|4.n(breast) ~ 0.03 cm"' 
and fi/ (breast) ~ 12 cm"'], see Heusmann et al., J. Biomed Opt. 1:425-34, 1996; Heusman 
et al., Proc. SPIE 2389:788-808. 1995; Kolzer et al., Proc. SPIE 2326:143-52, 1995. 

30 The experimental procedure to investigate the dependence of AI/Iq\^^^ on 

ji'joA}.i„ was based on the aixangement shown in FIG. 9. The laser diodes 10 were 
employed at the same six wavelengths used previously to determine the optical properties 
of the inclusion material and the background (752, 778, 786, 813, 830, and 840 rnn). In 
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5 this arrangement, the source fiber bundle 20 and detector fiber bundle 30 were arranged 
collmearly and then ends were fully immersed in the fomrerly characterized background 
medium to simulate an infinite geometry. The detector fiber bundle 30 fed the optical 
signals into the analyzer 40. The source-detector separation for these experiments 
remained a constant 6 cm, wliich is representative of the thiclaiess of a slightly 
1 0 compressed breast. Inclusions of different sizes and shapes were then suspended in the 
medium equidistant firom the source and detector using Pasteur pipettes to hold them in 
place. The pipettes were filled with background medimn to reduce the optical 
perturbation they might cause. 

The inclusion material was formed into in-egular shapes (by cutting cylindrical 
15 shapes with razor blades) and into a 10 cm-long cylinder using the two different 

mixtures, thus creating two different sets of the same sizes and shapes. The cylinder had 
a diameter of 1 .0 cm, while the irregular shapes were created to have the same volumes 
as spheres with 0.9 cm diameter (F= 0.43 cm^) and 1.4 cm diameter (F= 1.31 cm^), 
respectively. The source and detector optical fibers were scanned together over a 14 cm 
20 distance in the x-direction (as shown in FIG. 9) beguming at a distance of approximately 
7 cm fi-om the center of the inclusion in the x-direction. Scanning was performed at a 
rate of 0.65 mm/s and data was acquired every 0.88 s, providing a data point every 
572 \im- 

25 Methods 

Estimates of the optical signals were calculated using an analytical solution to the 
diffiision equation for a spherical inclusion embedded in an infmite turbid medium, see 

Boas et al, Proc. Natl. Acad. Sci. 51:4887-91, 1994. This solution was implemented in a 
PMI (photon migration imaging) software package developed by D. A. Boas et al. at 
30 Massachusetts General Hospital, Charlestown, MA, as part of pubhcly available 

software, http://www.imir.mgh.harvard.edu/DOT/toolbox.htm. Two types of theoretical 
analysis were performed. 
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5 The first analysis was aimed at deteimining whether the dependence of 

-^olmax °° l-^""' ^y^" is well approximated by a function of the product |j,jo'Ap.fl 

even outside of the perturbation regime, as in the cases of spheres having a diameter that 
is not much smaller than the source-detector separation and for values of A|j,„ tliat are not 
much smaller than ^„o. This assessed whether A/ /^olmax ' ^^''^ where Ap.^' = 0 

1 0 here, was independent of i^no and was not affected by |ijo' and A|a„ separately, but only by 
their product \isQ'L\ia. This analysis was carried out for two sphere diameters (1.4 and 
3.0 cm) in the case where the sphere is equidistant from the source and detector scanning 
lines (sphere center is 3.0 cm fi-om either fiber m the coUinear case). For the smaller 
sphere (1 .4 cm in diameter), a 1 .5-cm off-center position was also analyzed, where the 

15 distances firom the sphere center to tlie source and detector fibers were 1.5 and 4.5 cm 
respectively, in the coUinear condition. 

The second analysis was aimed at testmg the new methods for measuring tumor 
oxygenation, hi this analysis, the background absorption and scattering spectra, |a.„o(^.) 
and \x.sQ'm respectively, were set equal to typical absorption and reduced scattering 

20 spectra for breast tissue. The absorption spectra were determined by combining tlie 
absorption contributions firom typical concentrations of oxy-hemoglobin (-4.4 )j,M), 
deoxy-hemoglobin (-2.7 |xM), water (-39.5% v/v), and lipids (-60.5% v/v) in healthy 
breast tissue, as reported by Quaresima et al., Photochem. Photobiol. 62:4-14, 1998. hi 
the spectral region considered (680-880 nm), |a„o ranged from 0.023 cm"' (at 680 nm) to 

25 0.076 cm'' (at 880 nm). The background scattering spectrum was estimated fi-om data 
reported by Cubeddu et al, Appl. Phys. Lett. 74:874-76, 1999, for two healthy human 
subjects. In the spectral region considered (680-880 nm), ranged from 10.0 cm"' (at 
680 imn) to 7.8 cm"' (at 880 nm). 

In addition, spherical lesions with no scattering contrast (Ap,/ = 0) were also 

30 considered witli an absorption contrast provided by a hemoglobin concentration of 60 pM 
(corresponding to A|x„ = 0.12 cm"' at 800 nm). As in the first analysis, two sphere 
diameters (1.4 and 3.0 cm) in the centered case were analyzed (sphere half-way between 
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5 the source and detector scanning lines), along with one sphere diameter (1 .4 cm) in the 
off-center case (sphere 1.5 cm off-center). The maximum relative change in the optical 
signal induced by the spherical lesion ( A/ / /q ) was calculated at nine wavelengtlis 
(680, 720, 730, 758, 768, 776, 800, 840, aiid 880 mn) as a function of hemoglobin 
saturation values within tlie sphere over the range 0-100%. This choice of wavelengths 

10 covered the 680-880 nm diagnostic window and included the local maxima and minima 
in the deoxy-hemoglobin absorption spectrum (730 and 758 nm, respectively) and the 
local minimum in the oxy-hemoglobin absorption spectrum (680 nm). 

Testing of the Method 

15 FIG. 10 shows a table of the results of the measurements of the optical properties 

of the inclusions and background medium. The range of A|j,n is from 0.014 cm'' to 
0.12 cm'' when considering both sihcone mixtures, while the range of A|j,/ is jB-om -0.5 to 
-2.8 cm''. Given the background optical properties of the order of 0.02 cm"' for \i„o and 
9 cm'' for |iso', the absorption perturbations are about 70-600% of the background 

20 absorption, while the scattering perturbations, in absolute value, do not exceed 31% of 
the backgi'ound reduced scattering coefficient. These conditions are representative of the 
range of optical contrast offered by tumors in the human breast in vivo. 

FIG. 1 1 shows a typical result of the optical intensity (i) measured durmg a linear 
scan across the object location at x = 0. Because A|j,„ > 0, there is a decrease in detected 

25 intensity as the source-detector pair approaches the inclusion during the scan. The 

background value (/o) and tire maximum (in absolute value) intensity change lAI|max are 
indicated in FIG. 1 1 . The spectrum of (A/)max//o(?i'), is then used to guide the choice of 
the two wavelengths Xi and X2 that are used to measure the oxygenation of the embedded 
object according to Eq. (7). 

30 The experimental results for - AI ^hlmax ^ function of the product Hio'Ap,„ are 

reported in FIG. 12 for the irregularly shaped inclusions, and in FIG. 13 for the 
cylindrical inclusion. In both cases, the fact that the experimental data of M /-^olmax 
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5 a range of values of |j.jo' aiid Aj.i„ (see Table 1) are distributed along a single curve as a 
function of the product p,jo'A|J,„ is indicative of the dependence of AI/Iq\^^^ on f.i^o'A^n- 
In FIGs. 14 and 15, the functions of ii,o'A\x„ are represented by continuous lines. These 
experiments confinn that M //oUax °^ly ^ function of the product ^l,o'A^„ even for 
objects with a relatively large size, cylindrical or irregular shape, and relatively high 

10 absorption contrast. Furthennore, AI //oLax "^^^^ monotonically with H5o'Aja„. 

Figure 16 shows the calculated dependence of -AI/ /o|^^^ on the product |J,jo'A)i„ 
for a spherical object embedded in a uniform turbid medium. As described above, two 
spheres (diameters of 1.4 and 3.0 cm) were examined for the case where the sphere was 
equidistant from the source and detector fibers, and one sphere (diameter of 1 .4 cm) for 

1 5 the off-axis case where the sphere was 1 .5 cm off tlie mid-line between source and 
detector. The absorption coefficient of the background medium (\iao) is 0.06 cm"\ and 
the scattering perturbation (A^/) was set to zero. In all three cases considered, 
AI I /q Lax separately dependent on p,jo' and Ajj,„, but it only depended, 

monotonically, on their product (see FIG. 14). Furthermore, FIG. 15 shows that 

20 AJIIq was weakly dependent on n„o over the range of optical properties of interest. 

Computation of Oxygenation Levels 

The values of -All Iq\^^^ at nine wavelengths (680, 720, 730, 758, 768, 776, 
800, 840, and 880 nm) were computed for a spherical inclusion in a case that mimics a 

25 reaUstic condition in optical mammography. As described above, the background 
scattering and absorption spectra were set equal to representative spectra for healthy 
breast tissue obtained from literature data, see R. Cubbedu et al, Appl. Phys. Lett. 
74:874-76, 1999; Quaresima et al, Photochem. Photobiol. 67:4-14, 1998, wliich are 
incorporated by reference herein. Because of the results shown in FIG. 15, the 

30 background absorption spectrum plays a minor role in the determination of AT / /q | 



-27- 



wo 03/077750 



PCT/US03/07766 



5 The hemoglobin concentration of the embedded lesion was set to a value 60 Mgher 
than that in the backgi-ound, and the hemoglobin saturation of the sphere was varied over 
the range 0-100%. 

The resulting spectra of tlie product [iso'A\Xa for SO2 values of 0, 20, 40, 60, 80, 
and 100% ai-e reported in FIG. 16. On the basis of Eq. (3), the spectra of FIG. 16 are 

10 representative of the spectra of A/ //qL^ • The hemoglobin satui-ation in the different 
lesions were compared: (1) by using Eq. (5) and a ratio A^„(A.2)/A|.i„(?„) computed with 
perturbation theory [Eq. (2)] from the intensity change AI / measured at two fixed 
wavelengths; and (2) by using tlie new methods [Eq. (7)] at two wavelengths that are at 
least 40 mn apart and that minimize the absolute value of the difference 

15 A// /olmax '-'^2) " ^ ^ hlrj^^x ^^^^ ^ ' approaches take the background reduced 
scattering coefficient as input, which was Imown in our theoretical computations and 
which can be measured with time-resolved methods in a practical implementation to the 
human breast. The results of tlie perturbation analysis for a number of wavelength pairs 
are shown in FIG. 17 for a 1.4 cm-diameter sphere (on the ndd-line between source and 

20 detector), in FIG. 1 8 for a 3 .0 cm-diameter sphere (on the mid-line between source and 
detector), and in FIG. 19 for a 1.4 cm diameter sphere off the mid-line by 1.5 cm. The 
results of the new methods are shown in FIG. 20 (sphere diameter: 1 .4 cm; on the mid- 
luie), FIG. 21 (sphere diameter: 3.0 cm; on the mid-line), and FIG. 22 (sphere diameter: 
1 .4 cm; 1 .5 cm off the mid-line). 

25 The improved fit obtained to the actual oxygenation levels using the new methods 

is evident by comparing FIGs. 22-24 to FIGs. 19-21. None of the lines obtained using 
prior art perturbation analysis for a pair of wavelengths achieved a comparable fit to lines 
obtained using the new methods with a selected pair of wavelengths. These results show 
that the new methods provided accurate measurements of the oxygenation of spherical 

30 regions over tlie ftill range of oxygenation values and independent of the size and location 
of tlie sphere. The wavelength pairs that minimize the difference between A// 

"Imax 
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and A/ / hf^ll (which is the criterion used to select the particular wavelength pairs it 



in 



FIGs. 22, 23, and 24) for the cases of SO2 equal to 0, 20, 40, 60, 80, and 100% are 
indicated in FIG. 16. In FIG. 16, the wavelengtlis used are shown by the dashed lines, 
and the symbols identify the wavelength pairs that minimize the difference between 
A/ / /o |[^^ and A/ / Iq f^^ for each value of SO2. 

Additional Computational Metliods 

The new techniques discussed above quantify the oxygen saturation associated 
with A^fl , i.e., with the additional absorption at the tumor location with respect to the 
background tissue [see Eq. (5)]. These teclmiques are useful where (1) the absorption at 
the tumor location results from the sum of the background (healthy tissue) absorption 
{\i„o) plus the tumor contribution {[i/^ = A|j,„), or if (2) the tumor absorption is equal 
to |i„o + An« and A|in » |i„o (high-contrast tumor). In the case where |j,„^'^ = + A|ifl, 
but A\Xa is not much greater than |j.„o, the tumor saturation S02^'^ is more appropriately 
given by: 



where R = A\i„{X2)/A\i„(Xi). Because R is on Hie order of jiso'(^i)/Hso'(^2), it is typically 
close to 1. Consequently, S02^'^ given by Eq. (1 1) tends to the background saturation 
(S02^°^) in the limit A\i„ 0 and to the saturation based on A]ia (SO2) in the limit 
A|J,„ » |a„o. The derivative of S02^'^ with respect to A|j,„(X,i) at constant R is given by: 




[S Hb (^2 ) - e Hb02 (^ 2 )] + [S Hb02 (^ 1 ) - £ Hb (^1 )] ■ 




(11) 
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(dsol_] 

^^«0(^2)--^l-ta0(^l)][gHb(^2)£HbO2(^l)-gHh(^l)eHbO2(^2)1 

i^Eb (^2 ) - eHb02 (^2 )i^aO ) + ^l^a (^1 )] + [eHb02 )-^mih )]ko 0^2 ) + iJAu, (h )]f 

(12) 

Equation (12) shows that for any given value ofR, the sign of tlie derivative of 
SOz^'^ with respect to A\x„{h) is independent of Aia„(^i). As a result, S02^'^ monotonically 
increases (or decreases) from SOj^"^ to SO2 as the tumor conti-ast increases. Tlxis is 
gi-aphically shown in FIG. 23 for the case SOs^"^ = 70%, ^„o(800 nm) = 0.02 cm-' (a 
reasonable estimate of the background hemoglobin absorption in breast tissue), and SO2 
values of 0, 20, 40, 60, 80, and 100%. Therefore, the sign of S02^'^ - S02^°\ which 
indicates whether the tumor is more oxygenated or less oxygenated tlian tlie background 
healthy tissue, is accurately given by the sign of SO2 - S02^°^ independent of the tumor 
contrast. Furthermore, FIG. 23 indicates that for An„ > 0.1 cm'' and tumor saturation 
values greater than 40%, S02^'^ and SO2 differ by no more than -5%. Using an estimate 
of AyL„(k\), for instance on the basis of Eq. (2), can provide a refinement of the tumor 
saturation measurement from the values of SO2 (obtamed using the new methods) and 
S02^°^ (measured on healthy breast tissue). This use of such an estimate is appropriate 
where the tumor absorption is better represented by ]Xoo + A^a then by A|j,„. 

In Vivo Results 

The techniques and methods discussed above have been appUed to produce 
images of wj vivo tumors in human tissue. For example, FIGS. 24al to 24c2 show left 
cranio-caudal (Icc) and left oblique views (lob) of a 3.0 cm invasive ductal carcinoma in 
the left breast of a 53-year-old patient. The location of the lesion is known from x-ray 
mammography, and the nature of the lesion is known from biopsy or physical 
examination. FIGS. 24al and a2 show intensity (N) images of Ught transmitted through 
the breast. FIGS. 24bl and b2 show second-derivative (N") rniages of the transmitted 
image; and FIGS. 24cl and c2 show the relative oxygenation index of the features. 
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5 The second derivative image is generated by processing the N-image through a 

low-pass-filter (smoothing), calculating the spatial second derivative of N at each pixel in 
four directions (horizontal, vertical, two diagonals), talcing the minimum of the four 
directional second derivatives to enhance the detection of directional structures such as 
blood vessels, displaying on a gray scale tlie pixels associated with a negative second 

1 0 derivative (i.e., those corresponding to absorbance maxima), and setting tlie positive- 
second-derivative pixels to wMte. 

The cancer is detected in the left craniocaudal (Icc) view of the intensity image 
(FIG. 24al), together with a conspicuous blood vessel, but is not readily visible in the left 
obhque (lob) view (FIG. 24a2). The image conti-ast is enlianced in the second-derivative 

15 (N") images (FIGS. 24b 1 and b2), wliich display negative second-derivative values (gray 
areas) at the cancer location in both Icc and lob views. However, additional 
inhomogeneities (most likely associated with blood vessels) appear in the second- 
derivative images (24b 1 and b2) that were not visible in the intensity images (FIG. 24al 
and a2). The increased sensitivity afforded by the second-derivative images may be 

20 effectively complemented by the functional information provided by the oxygenation 

index images. Intliis case, the tumor location corresponds to areas of lowest oxygenation 
index in both the Icc and lob views (see FIGS. 24cl and c2). 

FIGS. 25a to 25c show left cranio-caudal (Icc) views of a 1 .5 benign mastopathy 
in the left breast of a 62-yeai--old patient. FIG. 25a shows an intensity (N) image of light 

25 transmitted through the breast. FIG. 25b shows a second-derivative (N") image of the 
transmitted image; and FIG. 25c shows the relative oxygenation index of the features. 
The mastopathy is visible in both the intensity image (FIG. 25a) and the second- 
derivative image (FIG. 25b). However, oxygenation level of the mastopathy is relatively 
high, as shown in the oxygenation level image (FIG. 25c), indicating that the tumor is 

30 benign. 

FIGS. 26a to 26c show a right obhque (rob) view of benign inflanmiation in a 
breast of a 65-year-old patient. FIG. 26a shows an intensity (N) image of light 
transmitted through the breast. FIG. 26b shows a second-derivative (N") image of the 



-31 - 



wo 03/077750 



PCT/US03/07766 



5 transmitted image; and FIG. 26c shows the relative oxygenation index of the features. 
The inflammation shows up as the dominant feature in the N-image as well as in the N"- 
image, but structural information is much more defined in the N" -image. Furthermore, 
the oxygenation index image shows the highest value of oxygenation in the region of the 
inflaiiamation, and lower oxygenation values elsewhere, indicating that the inflanimation 

10 is benign. 



OTHER EMBODIMENTS 
It is to be understood that while the uivention has been described in conjunction 
with the detailed description thereof, the foregoing description is intended to illustrate 
1 5 and not limit the scope of the invention, which is defmed by the scope of the appended 
claims. Other aspects, advantages, and modifications are within tlie scope of the 
following claims. For example, the new techniques can be used to image brain tissue and 
assess brain activity by measuring absolute levels of oxygenation, or to image muscle 
tissue and identify ischemic or underperfiised tissue regions. 
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5 WHAT IS CLAIMED IS: 

1. A system comprising: 

(a) a light source for illuminating a tissue sample; 

(b) a detector for detecting light transmitted tlirough the sample; and 

(c) a processor for calculating spatial second derivatives of products of sample 
10 thicknesses and intensities of the transmitted light. 

2. The system of claim 1, wherein the hglit source comprises a laser diode. 

3. The system of claim 1, wherein the hght source illuminates the sample with light 
havmg at least two wavelengths. 

4. The system of claim 1, wherein the light source illuminates the sample with liglit 
1 5 having at least one wavelength in tlie near infrared spectrum. 

5. The system of claim 1, wherein the detector comprises a photomultiplier tube or an 
avalanche photodiode. 

6. The system of claim 1, wherein the processor is a computer. 

7. The system of claim 1, wherein the spatial second derivatives, N'Xx^), are calculated 
20 from an equation N{x,y) = roaco/r{x,y)ac{x,y), wherein 

X andy are respective x- and y-coordinates of the coplanai* locations; 
ro is a maximum thickness of the sample; 

aco is the intensity at the location where the thiclaiess of the sample is ro; 
r{x^) is the thickness of the sample at location {x,yy, 
25 ac(x,y) is the intensity at location (x^); and 

N(x,y) is dimensionless number at a location (x,y). 
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5 8. The system of claim 1, further comprismg (d) a monitor to display aii image of the 
spatial second derivates. 

9. The system of claim 1, further comprising (d) a computer-readable medium to store 

values of the spatial second derivatives. 

10. A computer-readable medium comprising computer-readable instmctions for causing 
1 0 a computer system to : 

(a) receive data from a tissue sanrple that represent measures of tliicknesses of the 
sample and intensities of light transmitted through the sample; and 

(b) calculate spatial second derivatives of products of the thicknesses of tlie 
sample and intensities of the transmitted Ught to generate values that represent an image 

15 of the sample. 

11. The computer-readable medium of claim 10, further comprising computer-readable 
mstructions for causing the computer system calculate a spatial second derivative from an 
equation N(x^) = roaco/r{x^)ac{x,y) to generate values, N"(x^), wherein 

X and y are respective x- and y-coordinates of the values; 

20 I'o is a maximum tliickness of the sample; 

aco is the intensity at a location where the thickness of the sample is I'o; 
r(x,y) is the thickness of the sample at location (x^); and 
ac{x^>) is the intensity at location {x^f). 
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5 12. A computer-readable medium comprising computer-readable instructions for causing 

a computer system to: 



transmitted through the sample for at least two wavelengths of light; 

(b) select a pair of wavelengths to minimize a difference between relative changes 
1 0 in the intensities of the light ti-ansmitted tlirougli a region of the sample, wherein tlie 

relative changes in intensities of the Ught ai-e measured relative to a background intensity 
of light transmitted through the sample; and 

(c) process the relative changes in intensities of the Hght at tlie pair of 
wavelengths using measures of a background reduced scattering coefficient of the sample 

1 5 at the pair of wavelengtlis and molar extinction coefficients for deoxy-hemoglobin and 
oxy-hemoglobin to generate data values represent an absolute level of oxygenation of the 
region of tlie sample. 

13. The computer-readable medium of claim 12, further comprising computer-readable 
instructions for causing the computer system to generates values of the absolute level of 
20 oxygenation, SO2, of the region of the sample using a fonnula 



(a) receive data ixom a sample representing measures of intensities of hght 



SO2 = 



eHb(^^2)-eHb(^i)- 




Hb (^^2 ) ~ ^ Hb02 (^2 )3 + Hb02 (^1 ) ~ ^Hb (^1 )] ' 



wherein: 



25 



A/ / /fl 1^*''' is a maxunum value of a relative change in intensity at a first 
wavelength in the pair of wavelengths; 

A/ / /(, 1^^"^^' is a maximum value of a relative change in intensity at a second 
wavelength in the pair of wavelengths; 
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5 — p — ^ is a ratio of the background reduced scattering coefficient at the first 

wavelength to the background reduced scattering coefficient at the second wavelength; 

SHb(^i) and SHb(A.2) are the molar extinction coefficients for deoxy-hemoglobin at 
the first and second wavelengths; and 

subozih) and SHb02(^2) are the molai' extinction coefficients for oxy-hemoglobin 
10 at the first and second wavelengtlis. 

14. A method of detecting a region of a tissue sample having an optical absorption level 
different from a background optical absorption level, the method comprising: 

(a) obtaining thiclaiesses of the sample and intensities of light transmitted through 
1 5 the sample at a plurality of locations; 

(b) calculating spatial second derivatives of products of the sample thicknesses 
and the intensities of the transmitted light at the locations; and 

(c) comparing the spatial second derivatives to detect a region of the sample 
having an optical absorption level different from a background level of optical absorption 

20 in the sample. 

15. The method of claim 14, where the tissue sample is selected fi-om tlie group 
consisting of breast tissue, brain tissue, and muscle tissue. 

16. The method of claim 14, wherem the optical absorption is a characteristic of 
absorption of light by hemoglobin in tlie sample. 

25 17. The method of claim 14, wherein the region of the sample having an optical 

absorption level different from the background level of optical absorption within the 

sample comprises a lesion, tumor, or cancerous tissue. 

1 8 . The method of claim 14, wherein the transmitted light is generated by coplanar 
sources and wherein the plurality of locations are coplanar. 
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5 19. The method of claim 18, wherein the spatial second derivatives, N"{x,y), are 
calculated from an equation N{x,y) = roaco/r{x,y)ac{x,y), wherein 

X and y are respective x- and y-coordinates of the coplanar locations; 
ro is a maximum thiclaiess of the sample; 

aco is an intensity at the location where the tliickness of the sample is ro; 
10 J'(x,y) is a thickness of the sample at location (xj;); 

ac(x,y) is an intensity at location (x^); and 
N(x,y) is a dimensionless number at a location (x,y). 

20. The method of claim 14, further comprising: 
(d) illuminating the sample with light; and 

1 5 (e) detecting light transmitted through the sample at a plurality of locations. 

21 . The method of claim 20, wherein tlie hght comprises light having at least two 
wavelengths. 

22. The method of claim 20, wherein tlie light comprises light having at least one 
wavelength in the near infrared spectrum. 

20 23. The method of claim 14, further comprising: 

(f) displaying an image of the spatial second derivatives at the plurality of locations. 

24. The method of claim 14, further comprising: 

(f) storing values of the spatial second derivatives in a computer-readable medium. 

25 . A method of calculating an absolute oxygenation level of a region of a tissue sample, 
25 the method comprising: 

(a) selecting a pair of wavelengths of light to minimize a difference between 
relative changes in intensity of Ught transmitted through a region of the sample, wherein 
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relative changes in intensity of the Ught are measured relative to a background intensity 
of light transmitted thi-ough the sample; 

(b) obtaiiaing measures of background reduced scattering coefficients of the 
sample at the pair of wavelengths; and 

(c) calculatmg an absolute oxygenation level of the region of the sample using the 
relative changes in intensity of the Ught at the pair of wavelengths and the measures of 
the background reduced scattering coefficients of the sample at the pair of wavelengths. 

26. The method of claim 25, wherein the tissue sample is selected from the group 
consisting of breast tissue, brain tissue, and muscle tissue. 

27. The method of claim 25, wherein the region of the sample comprises a lesion, tumor, 
or cancerous tissue. 

28. The method of claim 25, wherein tlie pair of wavelengths of light are in the near 
infrared spectrum. 

29. The method of claim 25, wherein the difference bet\veen the relative changes in 
intensity is zero. 

30. The method of clami 25, wherein tlie absolute oxygenation level of the region is 
calculated using a formula 



^Hb 0^2)-~^Hh(\) 



^'.0(^2) M/I\[ 




[S Hb (^2 ) - £ Hb02 2 )] + [e Hb02 ) - £ Hb (^1 )] " 



wherein: 



SO2 is the oxygenation level of tlie region; 

A/ / /(, 1*^^^^^ is a maximum value of a relative change in intensity at a first 



-38- 



wo 03/077750 



PCT/US03/07766 



5 wavelength iii the p air of wavelengths ; 

^ ^ ^ » Imax ^ niaximmn value of a relative change in intensity at a second 
wavelength in the pair of wavelengths; 

^f^. 'v ^ °^ ^ background reduced scattering coefficient at the first 

wavelength to a background reduced scattering coefficient at the second wavelength; 
10 smih) and smih) aie molar extinction coefficients for deoxj'-hemoglobin at the 

first and second wavelengths; and 

EHboaC^w) and Emoiih) are molar extinction coefficients for oxy-henioglobin at 
the first and second wavelengths. 

3 1 . The method of claim 25, further comprising: 

1 5 (d) illuminating the sample with a plurality of wavelengths of light; and 

(e) detectmg light transmitted through tlie sample at a plurality of locations. 

32. The method of claim 25, further comprising (d) displaying an image of tlie absolute 
oxygenation level of the region at the plurality of locations. 

33. The method of claim 25, further comprising (d) storing values of tlie absolute 
20 oxygenation level of the region in a computer-readable medium. 

34. A method of deteimining whether a tumor in a tissue sample is malignant or benign, 
the method comprising: 

(a) selecting two wavelengths of light to minimize a difference between relative 
changes in intensity of light transmitted through the tumor for the two wavelengths, 

25 wherein the relative changes in intensity of the light for each wavelength are measured 
relative to a backgi-ound intensity of light transmitted tlirough the sample; 

(b) obtaining measures of background reduced scattering coefficients of the 
sample at the two wavelengths; 
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5 (c) calculating an absolute oxygenation level of the tumor in the sample by using 

the relative changes in intensity of the light for the two wavelengths and the measures of 
the backgi-ound reduced scattering coefficients of the sample for the two wavelengtlis; 
and 

(d) determining whether the tumor is malignant or benign based on the value of 
10 calculated absolute oxygenation level. 

35. The method of claim 34, wherein the tissue sample is selected j&om the group 
consisting of breast tissue, brain tissue, and muscle tissue. 

36. The method of claim 34, wherein the two wavelengths of light ai-e in the near 
infrared spectrum. 

15 37. The method of claim 34, wherein the difference between the relative changes in 
intensity is zero. 

38. The method of claim 34, wherein the absolute oxygenation level of the tumor is 
calculated using a formula 

1-1.0(^2) Ai/ir^ 

SO2 =— — -— , 

' ex AI 11 1 

[SHb(^2)-eHb02(^2)] + [eHb02(^,)-eHb(^,)]^^^f^ 

20 wherein: 

SO2 is the oxygenation level of the tumor; 

A/ / /(, j'^^^^ is a maximum value of the relative change in intensity at a first 
wavelength in the pair of wavelengths; 

A/ / /q fn^jx^ is a maxmium value of the relative change in intensity at a second 
25 wavelength in the pair of wavelengths; 
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— — — IS a ratio of the background reduced scattering coefficient at the first 

wavelength to the background reduced scattering coefficient at the second wavelengtli; 

f^Hhih) and SHb(X,2) are molar extinction coefficients for deoxy-hemoglobin at the 
first and second wavelengths; and 

£Hb02(^i) and SHb02(^2) are molar extinction coefficients for oxy-hemoglobin at 
the first and second w^avelengths. 

39. The method of claim 34, fuither comprising: 

(e) illuminating tlie sample with a plurality of wavelengths of light; and 

(f) detecting light transmitted through the sample at a plurality of locations. 

40. The method of claim 39, further comprising (g) displaying an image of the absolute 
oxygenation level at the plurality of locations within the tissue sample. 

41. The method of claim 34, further comprising, (e) storing values of the absolute 
oxygenation level of the tumor in a computer-readable medium. 

42. A method of determining whether a tumor in a tissue sample is maUgnant or benign, 
the method comprising: 

(a) obtaining thicknesses of the sample and intensities of light transmitted through 
the sample at a plurality of locations for two wavelengths of light; 

(b) calculating spatial second derivatives of products of the sample thicknesses 
and the intensities of the transmitted hglit at the locations for the two wavelengths of 
light; 

(c) calculating an oxygenation level of the tumor based on the spatial second 
derivatives for the two wavelengths of hght, the molar extinction coefficients of oxy- 
hemoglobin for the two the wavelengths of light, the molar extinction coefficients of 
hemoglobin for the two wavelengtlis of light, relative changes in intensity of the light for 
tlie two wavelengths of hght and the measures of the background reduced scattering 
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5 coefficients of the sample for the two wavelengths of light; and 

(d) calculating an oxygenation level of non-tumor regions of the tissue sample 
based on the spatial second derivatives for the two wavelengths of light, the molar 
extinction coefficients of oxy-hemoglobin for the two wavelengtlis of light, the molar 
extinction coefficients of hemoglobin for the two wavelengths of light, relative changes 

10 in intensity of the light for the two wavelengths of hght, and the measures of the 

backgi-ound reduced scattering coefficients of the sample for the two wavelengths of 
light; and 

(e) comparing the oxygenation level of the tumor with the oxygenation level of 
non-tumor regions of the sample to detennine whether the tumor is malignant or benign. 

15 43 . The method of claim 42, wherein the tissue sample is selected firom the group 
consisting of breast tissue, brain tissue, and muscle tissue. 

44. The method of claim 42, wherein the two wavelengths of light are in the near 

infrared spectrum. 

45. The method of claim 42, wherein the oxygenation level of a region of the sample is 
20 calculated using a formula 



0L = 



A[HbOJ- 



— , wherein 



A[Hb02]*+A[Hb]' 



A[Hb02]* = 
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A[Hb]*=. 







[i:4b02(^-)]-| 




[^ZsHb02(X,-)£Hb(^/)l 


1 




H 


^S£Hb02('^,)SHb(^,)j 


2 



and wherein: 

OL is the oxygenation level of the tumor; 

i is a wavelength index for the two wavelengtlis; 

SHb02 and snb are the molar extinction coefficients of oxy-hemoglobin and deoxy- 
hemoglobin, respectively; 

A[Hb02]* and A[Hb]* are relative values for the spatial changes in the concentrations 
of oxy-hemoglobin and deoxy-hemoglobin, respectively; and 

N" is a spatial second derivative of an intensity of transmitted Ught. 

46. The method of claim 42, fbrtlier comprising: 

(f) illuminating the sample with a plurality of wavelengths of light; and 

(g) detecting light transmitted tlirough the sample at a plurality of locations. 

47. The method of claim 42, further comprising, (f) displaying an image of oxygenation 
levels at the pluraUty of locations within the tissue sample. 

48. The method of claim 42, further comprising, (f) storing values of oxygenation levels 
at a plurality of locations within the tissue sample in a computer-readable medium. 
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FIG. 8 
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FIG. 9 
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FIG. 12 
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FIG. 17 
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FIG. 20 
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N" at 690 nm 




0 2 4 6 8 10 12 14 16 18 20 



FIG. 24B-1 




y(cm) 



FIG. 24B-2 



wo 03/077750 



18/20 



PCT/US03/077f.6 



Oxygenation Index 
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